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iN-^RODuc^ro:; 


As  indicated  by  the  title,  this  report  is  specifically  concerned  with  the 
lo: ''i tudinal  approach  dynamics  of  remotely-piloted  vehicles.  Certain  portions 
of  this  effort,  however,  are  felt  to  be  applicable  to  the  more  general  problem 
of  the  IFR  landing  of  aircraft  of  any  type.  In  this  report,  the  "guidance" 
function  specifically  refers  to  the  process  of  defining  the  longitudiral  tra- 
jectory of  the  aircraft  through  space  and  time  from  the  point  of  acguisition 
to  touchdown.  If  one  denotes  ■ to  be  the  angular  error  of  the  aircraft  rela- 
tive to  the  desired  approach  path  and  R the  horizontal  range  to  the  landing 
point,  then  the  desired  spatial  flight  path  may  be  described  by  • = f(R).  If 
R denotes  the  time  variation  of  the  horizontal  range,  the  trajectory  specifi- 
cation in  space  and  time  is. completed  by  R = g(R).  From  the  two  functional 
relationships  • = f(R)  a^^d  R = g(R),  one  may  uniguely  compute  the  force  vector 
■•eauired  to  fly  the  trajectory  through  space  and  time.  If  one  views  the  guid- 
ance function  as  specification  of  the  force  vector,  then  the  "control"  function 
is  ouite  clearly  the  process  of  developing  and  controlling  the  aerodynamic, 
oravi tational , and  inertial  forces  acting  on  the  airframe  to  equal  the  force 
sector  specified  by  the  guidance  function.  It  is  the  opinion  of  the  author 
tnat  the  quidance  formulation  contained  in  this  document  has  complete  applica- 
'ion  to  f.he  qeneral  problem  of  the  IFR  landing  of  aircraft  of  any  type;  how- 
ever. the  formulation  has  been  specialized  in  this  docum.^nt  to  the  case  of  a 
non-decelerati ng  (R  = constant)  approach  to  the  landing  point. 

The  "control"  function  is  strongly  influenced  by  the  nature  of  the  appli- 
cation. For  the  landinq  of  a pi  1 oted-ai rcraft,  the  control  function  emphasis 
will  invariably  center  upon  the  landing  display,  pilot  workload,  airframe  hand- 
linq  qualities  and  ot*^'''^  so-called  human  factors.  For  the  landing  of  a 
remotel v-pi 1 oted-vehi cl e , design  emphasis  will  tend  toward  the  auto-pilot  as 
ii.  is  tne  sole  means  tor  controlling  the  forces  on  the  airframe  in  the  pre- 
scribed manner.  Of  paramount  importance  to  these  two  superficially  different 
control  problems  is  the  specific  character  of  the  airframe  (e.g.,  helicopter, 
fixed  wing,  aspect  ratio,  drag  brakes,  wing  loading,  control  power,  density 
altitude,  etc.)  as  '.'ell  as  the  nature  of  the  on-board  sensors  (e.g.,  rate 
gyros,  accelerometers,  vertical  gyros,  etc.)  which  are  available  to  modify  the 
specific  character  of  the  airframe.  The  control  problem  considered  within  this 
document  is  the  automatic  landing  of  a mini-RPV  of  the  delta  wing  variety. 

The  guidance  and  control  of  the  lateral/directional  motions  of  the  RPV  are 
not  considered  here.  The  general  approach,  particularly  in  the  guidance  formu- 
lation area,  would  be  very  similar  to  that  used  for  control  of  the  longitudinal 
motions  of  the  aircraft. 

2.  GUIDANCE  FORMULATION 


The  genesis  of  any  guidance  concept  is  the  calculation  of  the  resultant 
forces  which  must  act  on  a vehicle  to  cause  it  to'  move  along  some  desired 
flight  path  at  some  speed  of  flight.  The  force  of  interest  here  is  the  one 
which  must  act  normal  to  the  aircraft  spatial  velocity  to  produce  the  desired 
flight  path  curvature.  Referring  to  Figure  1,  the  path  curvature  is  related 
to  the  required  incremental  acceleration  by 


a = 


V , 


V ^ 

dR  dt 


(2.1) 
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Figure  1 


r: 

It 


j 


Approach  coordinates 


F 


The  time  rate  of  change  of  ground  range  is, 

= R = cos  T (2.?) 

0 L. 

Note  tha*  R is  normally  a negative  number.  The  inclination  of  the  flight  path 
relative  to  the  horizontal  is  defined  by, 

tan  Y = (2.3) 

dR 


By  differentiating  (2.3)  with  respect  lo  the  ground  range  variable  (R),  one 
finds ; 


dR 


-cos2(-, ) 


d-  h 
dR2 


(2.4) 


Substitution  of  (2.2)  and  (2.4)  into  (2.1)  yields, 
a = v2cos2 i 

dD<^ 


(2.5) 


In  a typical  approach  maneuver  the  flight  path  angle  (y ) , the  desired  flight 
path  angle  (>p),  and  the  angular  error  (A)  are  small  enough  to  assume  SIN(  ) = 
TAN  ( ) = ( ) and  COS  ( ) = 1. 


(2.6) 


The  approach  guidance  problem  is  generally  more  suited  to  an  angular  equation 
formulation  rather  than  a linear  displacement  formulation.  The  conversion  may 
be  easily  accomplished  by  use  of  the  definition  of  the  angular  error  (X), 

> = - ,p  (2.7) 

R 


Differentiating  (2.7)  with  respect  to  the  range  variable  (R), 

dR  R dR  R 
and  then  once  again, 

o2>  ^ i _ 2 ^ + 2 /h  \ 
dR^  R dR2  r2  dR  R / 

Using  (2.3)  and  (2.7),  equations  (2.8)  and  (2.9)  become. 


(2.8) 


(2.9) 


(2.10) 
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And, 

1 ^n  2 

dR!^'  " R dR"'  ‘ R dR 

Rearranging  (2.11), 

d^h  _ R d-'  ^ 2 d1 
dR"  dR?  dR 

and  snhsti tution  of  (2.12)  into  (2.6)  yields, 

• 

a = R" 


dR‘ 


+ 2 


d 

dR 


1_ 


(?.n) 


(2.1?) 


(2.13) 


Expression  (2.13)  relates  the  acceleration  experienced  at  any  range  R.by  a body 
traveling  along  a shallow  flight  path  A = A(R),  to  its  ground  speed  (R)  and  to 
the  specific  nature  of  the  flight  path  (d^A/dR^,  dA/dR)  at  that  range  (R). 


Quite  often  one  encounters  a situation  where  a basically  responsive  ve- 
hicle must  first  be  controlled,  and  then  guided  with  precision  along  some  pre- 
determined flight  path.  The  "control"  of  the  vehicle  is  conventionally  associ- 
ated with  its  short  term  behavior  while  the  concept  of  "guidance"  relates  to  its 
more  long-term  behavior.  Quite  often  the  dynamics  of  the  control  of  the  ve- 
hicle are  substantially  faster  than  the  guidance  dynamics;  so  much  so,  that  the 
f i rst-ordered  analysis  of  the  two  functions  may  proceed  independently.  The 
landing  approach  of  an  RPV  is  of  this  class  of  problem.  One  may  take  advantage 
of  this  fact  by  assuming  a desirable  flight  path,  computing  the  guidance  com- 
mand required  to  follow  the  desirable  flight  path,  and  later  concern  oneself 
with  the  details  of  the  autopilot  required  to  follow  the  guidance  command 
(i.e.,  the  control  function). 

A stable  flight  path  A = A(R)  from  the  current  state  of  the  aircraft 
(A  = A.| , R = R-] ) to  the  desired  landing  point  (A  = 0,  R = 0)  should  satis  "y  the 
following  constraints: 


1.  When  R = 0,  A = 0,  and  dA/dR  = 0 


2.  When  R = Ri  , A = A , and  ^ 

dR 


.1  (Yi  + Ai  + yf) 


(2.14) 


Note  that  the  current  state  of  the  aircraft  is  designated  by  the  subscript  "1." 
As  four  arbitrary  coefficients  are  required  to  satisfy  the  four  endpoint  con- 
straints, assume  the  flight  path  to  be  of  the  well-behaved  form, 

n+1  n+2 

A = Co  + Cl , R + C2R  + Cg  R (2.15) 

The  quantity  n is  a constant  parameter  which  will  be  used  for  adjusting  system 
gain  at  acquisition.  The  constraints  at  R = 0 immediately  require  Cq  = C]  = 0. 
The  remaining  expressions  for  Cg  and  Cg  follow  from, 

'l  = C2  Ri""^^  + Cg  Ri"+2  (2.16) 
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■ i ('1  + -1  + ,r)  = (n  + 1)  C2  Rl"  + (n  + 2)  C3  r/  " (2.17) 

Simultaneous  solution  of  the  algebraic  equations  (2.16)  and  (2.17)  yields. 


C2  ■ R.|n  + 1 I (fi  3) 


^ '1  ^ 'FI 


(2.18) 


-(n  + 2)' 


1 'M 


(2.19) 


"he  nomi nal  flight  path  from  the  current  state  of  the  aircraft  ( R = Rl , 

• = , , = . .| ) to  the  landing  state  of  the  aircraft  (R  = 0,  >.  = 0,  , = ; p)  is. 


/n  \n  + 1 - 


(n  + 3)  ']  + i] 


+ . r + 


n + 2 


|^-(n  + 2)  - Yi  - ':p1 


(2.20) 


Also, 


= (n  + 1) 


R 1 \ (n  + 3)  A,  + Yi  + vp 

W,l  [r,)  ^ J 


(n  + 2) 


R + 1 /I 


(n  + 2)  - Yp! 


(2.21) 


iV  ’ ^ l(n  + 3)  + Y-,  + yF^  + 


' n / 

(n  * 1)  (n  + 2)(^J(i)  A, 


'1  ■ 3p  , 


Substituting  (2.21)  and  (2.22)  into  (2.13)  yields. 


a = r2 


|(n  + 1)  (n  + 2)^  jjn  + 3)  A.,  + y-,  + Yf)j 


(2.22) 


(2.23) 


(n  + 2)  (n  + 3)/  ^ ^ i_  ^(n  + 2)  A.,  - Y^  - Y 


5 


Again,  (2.23)  describes  the  acceleration  profile  required  to  fly  the  nominal 
flight  path  from  the  current  state  (R  = R-i,  , = y-|  , A = >■))  to  the  landing 
state.  The  current  coniiiand  (ac)  results  from  using  R = R],  in  (2.23).  For 
current  guidance  command  generation  purposes ,• the  subscripts  "1"  are  dropped 
as  all  variables  are  current.  The  resultant  guidance  command  is, 

ac  (n  r 2)  -(n  + 3)  - 2,  - 2^p  ‘ (2.24) 

R - 


In  effect,  *:hp  aircraft  never  attempts  to  return  to  any  previous  nominal  flight 
path,  but  instead,  continually  recomputes  a new_nominal  path  from  the_current 
state.  As  the  measurable  quantities  are  ■ and  A , use  of  (2.10)  and  A = R d.-/ 
dR  yields. 


A 


-R 


( ‘ + A + I p ) 


Use  of  (2.25)  in  (2.^,4)  yields. 


(2.25) 


ac  = (-  2R)  (n  + 2)  (A^  - A) 


where 


(2.26) 


(2.27) 


Expression  2.27  will  be  used  in  a later  section  to  command  limit  A and  hence 
the  maximum  descent  angle  (y).  The  acceleration  command  of  (2.26)  can  be  ma- 
nipulated into  a second-ordered,  linear  differential  equation  by  use  of  (2.13). 


r2  ( R ^ + 2 d_^\,  _2R  (n  + 2)  TRJiLjLJl  _ - R (2.28) 

\ ^ dR/  L 2 R dR.J 

Rearrangement  of  (2.28)  yields, 

r2  ^ . 2 (n  + 1 ) R -^  + (n  + 1 ) (n  + 2)  A = 0 (2.29) 

rip?  dR 


The  particular  solution  to  2.29  which  satisfies  the  desired  boundary  conditions 
is  obviously  equation  2.20.  (If  this  is  not  obvious,  the  reader  should  recall 
that  the  solution  2.20  was  the  assumed  starting  point. ) Equation  2.29  will  be 
used  in  the  next  section  on  limiting  system  gain. 

Figure  2 depicts  the  response  of  2.26  for  acquisition  of  the  glide  slope 
(-,'p  = 40)  at  R = 3,000  feet,  with  an  initial  altitude  error  of  100  feet  (A  = 
1/30  radian).  The  range  and  altitude  error  channels  are  dynamically  rescaled 
at  R = 300  feet  for  accuracy.  The  inherent  stability  of  2.26  for  small  values 
of  R is  quite  clear. 
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3.  SYSTEM  GAIN  CONTROL 

The  quidance  concept  of  section  2 assures  that  > = / = 0 as  R goes  to 
zero,  by  increasing  the  system  gain  to  infinity  as  R goes  to  zero.  Ideal  sys- 
tem stability  is  insured  by  the  functional  relationship  of  > = [>(R)]  to  R as  R 
goes  to  zero  (i.e.,  LIMr  .q('^/R)  = 0.  Introduction  of  typical  dynamic  lags  in 
aircraft  response  (i.e.,  a ^ a^)  results  in  stability  problems  somewhere  near 
the  landing  point.  The  smaller  the  response  lags,  the  nearer  the  aircraft  can 
approach  the  landing  point,  before  stability  problems  result.  The  intent  of 
this  section  is  to  modify  2.26  by  incorporating  a near-field  gain  control.  By 
adjustment  of  the  gain  limit,  the  guidance  law  may  be  made  compatible  with  the 
dynamics  of  any  aircraft/autopi lot  control  system.  As  usual,  the  price  one 
pays  for  adequate  system  stability  is  degraded  performance. 

The  author's  initial  attempt  at  gain  control  was  to  use  an  altered  form  of 
2.29  to  control  the  aircraft  in  near  field  (R  ^ R^).  That  is, 

R^  - 2 (n  + 1)  Rm  ^ + (n  + 1)  (n  +2)  A = 0 (3.1) 

m dR 

The  above  equation  is  linear,  with  constant  coefficients.  As  such,  the  most 
obvious  source  of  increasing  gain  has  been  eliminated.  Using  2.13  in  3.1  one 
fi nds ; 

a^  = 2 (n  + 2)  R > - (n  + 1 ) (n  + 2)  R^  (X/R)  (3.2) 

^ - ^m 

ac  = 2R  X+/lV2  (n  + 1)  R X -(n  + 1)  (n  + 2)  R^  x (3.3) 

\^m/U  Rm  J 

Equation  3.2,  which  represents  the  far-field  guidance  law,  is  identical  to 
2.26.  Equations  3.2  and  3.3  are  identical  for  R = R^  so  a smooth  transition  is 
assured.  The  near-field  law  (3.3)  is  precisely  equivalent  to  the  differential 
equation  (3.1).  Its  odd  form  reflects  the  fact  that  the  autopilot  will  control 
normal  acceleration  (d^h/dR^)  rather  than  (d^X/dR^). 

The  aircraft/autopilot  system  contributes  a relatively  constant  dynamic 
lag  in  d^h/dR^.  As  apparent  from  2.11,  the  importance  of  this  lag  to  the  ve- 
havior  of d^x/jR^ , increases  with  decreasing  range.  Although  3.1  removes  the 
most  obvious  source  of  increasing  gain  in  the  near  field,  the  fact  that  the 
autopilot  can  only  control  d2h/dR2  and  not  d^X/dR^,  reouired  another  approach. 

Returning  to  2.29,  consider  a transformation  of  variables  defined  by 


Substitution  of  3.5  and  3.6  into  2.29  results  in. 


o d^h 


2 (n  + 2)  R ^ + (n  + 2)  (n  + 3)  = 0 

dR 


(3.7) 


''he  above  equation  describes  the  ideal  response  of  the  guidance  law  in  terms  of 
altitude  error  relative  to  the  desired  straight  line  approach  path  of  inclina- 
tion ,p.  It  is  completely  equivalent  to  2.29  which  relates  system  response  in 
terms  of  the  angular  error  (A)  rather  than  the  altitude  error  (h^-). 

Near  field  guidance  will  be  based  upon  the  following  linear  differential  equa- 
tion with  constant  coefficients  obtained  from  3.7, 


- 2 (n  + 2)  Rm  ^ + (n  + 2)  (n  + 3)  h = 0 
m 2 dR 


(3.8) 


Since  d^h/dR^  is  equal  to  d^h^/dR^,  the  first  term  is  measurable  and,  in  fact, 
under  control  of  the  autopilot.  The  near-field  guidance  law  formulation  will 
be  completed  by  elimination  of  dh^/dR  and  hg  in  favor  of  the  measurable  quanti- 
ties ( and  A. 


Using  3.4  and  3.5  in  3.8,  ore  finds. 


(3.9) 


R .<  Rm 


ar  = - 2 (n  + 2)  (R)/ L-Va.  - 


(3. in) 


^1  - (n  ^ _3)Y L 
i 2 \ Rm 


(3.11) 
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Expressions  3.10  and  3.11  are  equal  to  the  far-field  guidance  expressions  2.26 
and  2.27  when  R = Rm  insuring  a smooth  transition.  A single  guidance  expres- 
sion encompassing  both  the  near-field  and  far-field  laws  is. 


ac  = (-  2R)  (n  + 2)/  i_\  (X^ 


- M 


and. 
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^ / 

c = (-  R) 

1'  -1 

V 2 j\' 

R"  = R 

if 

R >R^ 

R = Rm 

if 

R R^ 

(VR) 


(3.12) 


(3.13) 


The  guidance  law  operates  in  the  far-field  on  the  measured  angular  information 
with  increasing  gain  as  the  range  decreases.  At  some  range  (R  = R^)  the  guid- 
ance concept,  in  effect,  converts  the  measured  angular  information  (X,  X)  to 
its  linear  equivalent  (he,  ii^)  and  then  freezes  the  system  gain  at  the  R^ 
value.  If  the  system  is  stable  at  a range  R^,,  it  will  be  stable  for  all  R $ Rm- 
One  should  realize,  however,  that  the  system,  will  respond  to  a gust  input  which 
occurs  at  any  range  R < Rm  with  the  same  urgency  as  if  it  had  occurred  at  R = 

Rm.  One's  obvious  goal  then,  is  to  make  sure  that  the  system  achieves  its  max- 

imum tolerable  gain  at  R = Rm  to  minimize  the  effects  of  near-field  disturbances 
due  to  wind  gusts  and  manual  tracking  errors. 

To  quantify  the  near-field  gust  response  of  this  guidance  formulation, 
consider  the  response  of  3.8  to  an  aerodynamical ly-produced  disturbance.  Later 
in  this  document,  the  effects  of  the  autopilot  design  upon  gust  alleviation/ 
amplification  will  be  addressed.  There,  one  will  find  that  the  pitch  attitude 
behavior  of  the  stabilized  airframe  in  response  to  a normal  wind  gust  (AWg) 
disturbance  controls  the  degree  of  gust  amplification/alleviation  (i.e.,  the 
resulting  h or  y response).  If  the  airframe  does  not  pitch,  a wind  gust  of 
magnitude  AWg  will  soon  translate  the  airframe  at  that  rate  (i.e.,  li  = AWg). 

The  time  lag  to  achieve  this  rate  is  predominately  the  "short  period  mode"  of 

the  airframe.  For  the  analysis  to  follow,  this  time  lag  is  ignored  and  equat- 
ion 3.8  is  subjected  to  the  following  initial  conditions  when  t = Oj 

he  = 0,  = AWg  (3.14) 

dt 

where  AWg  is  positive  for  an  up-gust.  The  following  change  of  variables  is 
convenient  for  analytical  purposes. 

P>T  - Rq  " R (3.15) 

where  Rp  is  the  range  from  the  desired  landing  point  at  which  the  wind  gust 
occurs,  and  Rj  is  the  transposed  range  variable.  Equation  3.8  becomes. 
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^ + 2 (n  + 2)  Rm  ^ + (n  + 2)  (n  + 3)  he  = 0 
dRj2  dRj 

subject  to  the  initial  conditions,  at  Rj  = 0, 

he  = 0,  :/  AWg ^ 

dRT  I--) 


(3.16) 


(3.17) 


The  variable  Aag  is  physically  the  change  in  airframe  angle  of  attack  produced 
by  the  wind  gust  (i.e.,  ^'Y(O)  = Aag  due  to  A9  being  assumed  equal  to  zero). 
Since  equation  3.6  is  linear  with  constant  coefficients,  its  solution  is 
easily  obtained  by  Laplace  transform  techniques.  Using, 


^rf^(t)l  = F(s)  - S^"''  f(0)  - f"(0).... 

The  Laplace  transform  of  the  solution  is 

A^g 

^^e(S)  = 5^  + 2 (n+2)  S + (n  + 2)  (n  + 3) 


(3.18) 


(3.19) 


The  nondimensional  solution  is. 


(3.20) 


(3.21) 


Figure  3 is  a plot  of  (3.20)  for  a family  of  values  of  the  parameter  (n).  This 
figure  indicates  that  with  Rm  held  constant,  increases  in  the  parameter  (n)  de- 
crease the  maximum  error  amplitude  (hf).  The  worst  case  location  Xo*"  'the  ar- 
resting gear  (R  = 0)  is  where  the  peak  gust  response  occurs  C^j  = Rjmax)  any 
value  of  the  parameter  (n).  As  indicated  previously,  the  intent  is  to  cause  the 
system  gain  to  be  at  its  maximum  tolerable  value  at  R = Rm.  This  in  effect 
specifies  a single  value  for  n so  that  conclusions  made  concerning  variations 
in  n,  with  Rm  held  constant,  are  invalid.  A change  in  n demands  a change  in  Rm 
to  hold  the  system  gain  constant.  Consider  the  undamped  natural  frequency  of 
oscillation  (l«l^)  and  damping  factor  (E^)  of  3.16. 


W = / (n  + 2)  (n  + 3) 

^ R 

^m 

Eh  = /rTt-  2 
/n  + 3 


(3.22) 

(3.23) 


n 


1 


If  one  makes  an  a’priori  assumption  that  R % lOOn  feet  and  n lO  are  typical, 
then  one  may  approximate  (3.22)  and  (3.23)  by. 


(3.24) 


(3.25) 


The  frequency  of  oscillation  (Wh)  best  characterizes  the  concept  of  system  gain 
and  stability  so  that  the  desire  is  to  examine  the  behavior  of  (3.20)  with 
held  constant  rather  than  R^. 

The  peak  value  for  hr  (Rj)  results  from  first  differentiating  (3.20)  with 
respect  to  R-[-,  setting  the  result  equal  to  zero,  and  solving  the  resultant 
equation  for  the  value  of  Rj  at  which  the  peak  occurs  (RjmAx) 


Rtma:.  = ^ ^ 

I n + 2 


(3.26) 


Again  assuming  n ~ 10, 


Rtmax  = L_ 

n + 2 


(3.27) 


K,-  = 

+ 2 


(3.28) 


Substituting  (3.27)  into  (3.20)  one  finds  the  peak  nondimensional  amplitude 
(*^EMAX^ 


^EMAX  " p-'’ 

V n + 2 SIN 


(3.29) 


Again  using  n '\.  10 

h - e-1 

TMAX  ^ 

n + 2 


1 ,,.-1 


hgn^v  = RimO  = e~  Wh 
.’-x_  n + 2 


(3.30) 
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For  the  autopilot  control  system  discussed  in  a later  section,  stability  con- 
siderations dictate  that  the  minimum  value  of  Wh"^  is  on  the  order  of  60  feet. 
From  (3.30),  the  worst  case  gust  response  for  the  RPV  system  of  this  document 
is  app>'oximately  21  feet  of  error  per  radian  of  Aag  or  at  an  approach  speed  of 
85  feet  per  second,  approximately  1/4  foot  of  error  per  foot  per  second  of  AWq. 


As  indicated  by  (3,28)  and  (3.30),  the  gust  characteristics  of  the  guidance 
law  during  the  final  approach  phase  (R  ^ R^)  are  only  a function  of  the  chosen 
spatial  Period  (W^"^)  of  the  system.  It  is  suggested  here  that  the  choice  of 
the  particular  combination  of  n and  R^  to  achieve  this  chosen  spatial  period 
(Wh"'  0.  60  feet)  be  influenced  by  two  considerations.  These  are  (1)  the  amount 
of  time  (-  Rni/ft)  desired  for  high  gain,  concentrated  tracking  in  the  final  ap- 
proach phase  and  (2)  the  system  gain  at  acquisition  (R  '■  3000  feet).  For  the 
present,  judgment  (Figure  4)  was  used  to  select  the  Rm  = 1000  feet,  n = 15  op- 
tion for  this  brief  automated  simulation  test.  This  allows  about  12  seconds 
for  the  final  approach  phase  and  reduces  acquisition  system  again  (feet  per 
second^  of  commanded  acceleration  per  radian  of  measured  angular  error  \)  by  a 
factor  of  3 from  its  value  at  initiation  of  the  final  approach  phase. 


Figures  5 through  8 depict  the  gust  response  of_3.13  without  gain  limiting 
(i.e.,  Rpi  = 0),  when  subjected  to  an  upgust  (AWg  = -RAaG  = 5 FPS)  at  several 
ranges.  The  system  is  initially  trimmed  on  the  descent  path  and  is  perturbed 
from  this  path  by  the  upgust.  Since  the  assumption  of  A9  = 0 and  initially 
h = AWq  are  retained,  the  associated  step  change  in  * is  clear.  Since  the 
system  gain  is  continual ly  increasing,  the  peak  Ah^  decreases  for  decreasing 
range.  A1  thougiTTTgure  3 is  only  truly  valid  for  values  of  R less  than  R^,  use 
of  R^,  = Rd  gives  a good  approximation  to  the  results  of  Figures  5 through  8 
where  R^  - 0.  For  example,  from  Figure  5,  the  peak  altitude  error  is  Ahc  = 4 
feet.  From  Figure  3 for  Rm  = Rq  ~ feet,  the  predicted  value  is  Ah£  = 3.7 

feet  (Ah£  = 0.021  Auq  RO).  As  indicated  by  the  nondimensional  results,  the  peak 
h£  remains  approximately  proportional  to  the  range  (Rp)  at  which  the  excitation 
occurred  (e.g.,  4 feet  at  Rq  = 3000  feet,  3 feet  at  Rq  = 2250  feet,  2 feet  at 
R[)  = 1500  feet  and  1 foot  at  Rq  = 750  feet  as  shown  by  Figures  5 through  8). 

Figures  9 through  12  depict  the  ideal  system  gust  response  with  gain  limit- 
ing. Figures  9 and  10  were  computed  with  Rm  = 3000  feet  and  demonstrate  the 
desired  invarient  Ahf  response  for  all  values  of  R;;  Rm  . These  two  traces  should 
be  compared  with  Figure  5 (Rm  = 0)  where  the  disturbance  was  applied  at  Rp  = 

3000  feet.  Figures  11  and  12  were  computed  with  R^,  = 1000  feet  and  demonstrate 
the  ideal  system  response  for  the  suggested  n = 14,  Rm  = 1000  foot  design  con- 
figuration. For  the  gust  response  of  the  design  configuration  for  R '■  R„^,  one 
should  refer  to  Figures  5,  6,  and  7 as  gain  limiting  is  not  in  effect  for  R 

Again,  it  should  be  emphasized  that  the  preceding  results  for  A9  = 0 repre- 
sent the  worst  case  as  far  as  gust  sensitivity  is  concerned.  When  the  airframe 
is  free  to  respond  in  pitch  attitude  to  a AWq  disturbance,  the  peak  Ah£  will 
decrease. 


-I 
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Figure  4.  Spatial  period  of  normal  acceleration  guidance  law 


15 


Figure  8.  Ideal  system  gust  response  (AW, 
Rr,  = 750  ft,  n = 15) 
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Figure  9.  Ideal  system  gust  response  (AWg 
3000  ft,  Rd  = 2250  ft,  n = 15) 
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4.  PATH  ANGLE  CONTROL 


by. 


The  power  required  to  sustain  the  flight  of  an  aircraft  may  be  represented 


= L • W.  + D • Vw  + Vw  • m ^ • mg 

dt 


(4.1 


The  first  term  is  the  induced  power  required  to  sustain  lift  (L  ' Wj).  The 
second  term  is  the  power  required_to  overcome  the  aerodynamic  drag  of  the  air- 
craft at  the  airspeed  of  flight  (D  ' \l^).  The  third  term  is  th^e  power_re- 
quired  to  inertial ly  accelerate  the  aircraft  in  any  direction  (V  • m dV/dt) 
such  as  turning,  pull-up,  or  changing  the  speed  of  flight.  The_iast  J^erm  is 
the  power  required  to  sustai n the  orientation  of  the  airspeed  * mg)  rela- 
tive to  the  gravitational  force  (e.g.,  a steady  climb).  The  source  of  the 
power  input  to  the  system  is  the  engine  and  it  is  assumed  Preo  = PenG  (''•e., 
the  engine  instantaneously  responds  to  the  flight  power  demanas).  The  minimum 
value  for  the  righthand  side  of  4.1  is,  therefore,  the  flight  idle  power  output 
of  the  engine. 

When  the  aircraft  acquires  the  glide  slope,  the  autopilot  is  able  to  pitch 
the  aircraft  over  to  begin  ;the  descent  and  also  decrease  Pe[^g  main_tain  the 
desired  approach  airspeed  (Vw).  As  the  orientation  of  the  airspeed  {M^)  steep- 
ens  (V^^_»mg  mo_re  negative),  less  and  less  power  is  required  to  maintain  air- 
speed *m  dV/dt  H 0).  When  fli_ght  i^le  power  is  reached,  further  increases 
in  the  steepness  of  the  airspeed  (\^  • mg),  must  be  accompanied  by  an  increase 
in  the  magnitude  of  th£  airspeed  (V^^).  The  cleaner  the  aircraft  is  in  the  ap- 
proach configuration  (D  • Vw),  the  smaller  will  be  this  limiting  path  angle. 
These  considerations  lead  to  the  conclusion  that  the  guidance  law  must  incor- 
porate a means  for  positive  control  of  the  aircraft  flight  path  angle  (y)  in  a 
descent.  As  previously  indicated,  this  is  necessary  to  avoid  undesirable  in- 
creases in  the  approach  airspeed  which  must  be  arrested  at  touchdown. 

To  control  path  angle  from  either  the  air  or  the  ground,  one  must  be  able 
to  directly  measure  it  or  to  be  able  to  at  least  compute  it  from  other  measure- 
ments. Since  the  only  guidance  measurements  are  X,  X,  and  R,  the  latter  ap- 
proach is  required.  Recalling  2.25, 


X = R ^ = zi  (^  + + X) 

dR  R 


(4.2) 


The  only  unknown  quantity  in  4.2  is  the  desired  flight  path  angle  (y).  If  the 
desired  maximum  descent  angle  (-  Y|  ) is  substituted  into  4.2,  then  a limit  on 
the  allowable  \ results. 


'Vlim 


- R (-  Y + Y,  + X) 
R 


(4.3) 
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The  intent  is  to  subject  3J2_and  3.13  to  a command  limiting  process  as  defined 
by  4.3.  The  command  limit  (-^clIM)  is  variable  being  a function  of  instantaneous 
range  and  angular  error.  In  addition  to  functioning  as  a command  limit.  An  im 
can  also  act  as  an  automatic  abort  signal.  Normally  is  negative  as.  'j  is 

certainly  qreater  than  T'p  and  A is  nominally  zero.  A positive  value  for  aClIM 
indicates  that  >.  has  become,  either  at  acquisition  or  later  in  the  approach,  too 
large  to  execute  the  maneuvers  to  touchdown.  This  would  be  even  more  important 
if  a positive  limit  was  imposed  on  the  guidance  law  to  avoid  attempting 

to  climb  at  too  steep  an  angle  for  the  maximum  engine  power  available.  The  pos- 
itive limit  will  not  be  considered  further  here  as  a complete  engine  simulation 
is  necessary  to  explore  this  boundary  of  the  flight  envelope. 


The  complete  guidance  law  is. 


4c 


(-  2r)  ■) 


and, 

R ' = R i f R >.  Rm  1 

R = Rm  i f R < Rm  [ 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


Figure  13  depicts  the  ideal  system  response  (Aa  = ) for  the  guidance 

law  (4.4)  from  acquisition  at  R = 3000  feet  and  Ah^  = 100  feet  to  touchdown. 

The  channel  labeled  "LOGIC'.'  indicated  the  fraction  of  the  flight  path  during 
which  command  limiting  of  was  in  effect  (>j  = 8^). 

5.  AIRFRAME/AUTOPILOT  REPRESENTATION 

There  is  a tendency  to  consider  the  complete  airborne  computing  capability 
(analog  or  digital  considerations  are  of  no  consequence)  as  the  "auto-pilot." 
Linder  such  circumstances,  the  design  objectives  of  the  guidance  function  and 
the  control  function  may  become  unnecessarily  interrelated  or  even  obscured. 

The  view  taken  within  this  report  is  that  the  control  function  pertains  to  the 
control  of  forces,  of  any  origin,  which  act  on  the  aircraft.  The  aerodynamic 
and  inertial  forces  on  the  aircraft  are  predominately  a function  of  the  trans- 
lational (U,  V,  W)  and  rotational  (p,  q,  r)  rates  of  the  airframe  and  their 


derivatives,  while  the  gravitational  force  is  a function  of  the  angular  orien- 
tation (({>,  P)  of  the  airframe.  Within  this  framework,  control  of  airspeed, 
altitude  rate,  glide  path  angle,  pitch  and  roll  attitude,  and  normal  accelera- 
tion are  all  considered  autopilot  functions,  while  control  of  heading  {\p) , 
altitude  (h),  translational  position  (N,  E)  or  even  the  ’ variables  of  this 
document  are  all  considered  guidance  functions. 

The  pertinent  quantities  for  autopilot  control  in  this  document  are  pri- 
marily the  normal  acceleration  and  secondarily,  the  airspeed.  The  design  ob- 
jectives for  the  autopilot  are: 

a.  To  stabilize  the  short  period  and  phugoid  modes  of  airframe  motion. 

b.  To  minimize  the  variability  in  normal  acceleration  response  per  unit 
of  command  input. 

c.  To  minimize  the  change  in  flight  path  angle  (y)  produced  by  a vertical 
wind  gust. 

The  goal  of  this  section  is  the  synthesis  of  an  autopilot  concept  rather 
than  an  autopilot  design.  No  attempt  is  made  to  optimize  the  parameters  with- 
in the  concept  to  best  satisfy  the  above  objectives.  Emphasis  is  placed  upon 
qualitatively  showing  how  each  parameter  impacts  the  various  objectives  and 
upon  arriving  at  a conservative  design  framework  for  later  optimization.  The 
technical  approach  used  herein  relies  heavily  upon  analog-computer  simulation 
as  a design  tool.  Analysis  is  used  to  explain  the  results  of  the  simulation. 

The  simulation  of  this  report  utilizes  a linearization  of  the  airframe 
dynamics  about  a nominal  condition  of  straight  and  level  flight  at  the  desired 
approach  airspeed.  The  stability  derivatives  are  estimates  of  a typical  flying 
wing  RPV  configuration.  The  specific  characteristics  of  the  engine  have  been 
neglected  for  this  effort  and  only  the  response  of  the  airframe  to  a perturba- 
tion in  thrust  (AT/m)  is  considered.  While  perturbation  aerodynamics  are  used 
throughout,  most  kinematic  relationships  were  retained  in  their  nonlinear  form. 
The  coordinate  system  used  herein  is  of  the  body  fixed  variety  with  its  origin 
at  the  center  of  mass  of  the  airframe.  The  linear  accelerations  along  the  ro- 
tating X axis  (i'  ) and  Z axis  (w) , and  the  angular  accelerations  about  the  Y axis 
(0)  are; 


AU  = Xw  AW;,  + Xu  AU  - g A0  + AT/m 

(5.1) 

• • • • 

AW  = U A0  - Aa  = U (A0  - Ay) 

(5.2) 

A6  = Mg  A0  + AWa  + M5  A6^ 

(5.3) 

where 

Aa  = - Z|5  A6f;  - AW - Zu  AU 

(5.4) 

AWa  = AW  + awg 

(5.5) 

U = Up  + AU 

(5.6) 
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The  quantity  Aa  is  the  perturbation  normal  acceleration  of  the  center  of 
mass.  A normal  accelerometer  mounted  a distance  ifi  front  of  the  airframe 
center  of  mass  will  sense. 


Aa^  = Aa  + 


A0 


(5.7) 


The  guidance  law  functions  in  a closed-loop  fashion  on  the  angular  error  (A) 
hence  from  (2.25), 


or 


A - [Y  + Yp  + A): 


A = -i  [ - AW  + U (A0  + Yp  + A)  ] 
R 


and 


A = 


A dt  + A, 


(5.8) 


(5.9) 


(5.10) 


The  altitude  error  at  any  range  is  found  from  the  open- loop  calculation 

hp  = R A (5.11 ) 

so  that  one  must  be  especially  cautious  of  scaling  difficulties  (e.g.,  A and 
R are  both  app>^oaching  zero).  The  simulation  results  within  this  document  may 
be  suspect  for  perhaps  the  last  2 seconds  of  flight  (170  feet)  due  to  analog 
scaling  difficulties  with  A and  R.  This  is  purely  a simulation  difficulty  and, 
in  any  event,  determination  of  absolute  errors  at  recovery  is  not  the  intent 
of  this  effort. 

The  following  stability  derivatives  were  used  in  the  simulation: 

X„  = .2  (Sec  -■') 

Xjj  = -.1  (Sec  -'') 

= -3.  (Sec  "1 ) 

Zu  = -.75  (Sec  ) (5.12) 

Z.j  = -50.  (Ft  - Sec  ~^) 

Mw  = -.1  (Ft  - Sec  ■'') 

Me  = -.6  (Sec  ■'') 

= -20. (Sec  -2) 

Ug  = 85.  (Ft  - Sec  "I ) 


Appendix  A contains  the  scaled  analog  computer  simulation  diagrams  for  equat- 
ions (5.1)  through  (5.11). 


k 
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The  set  of  equations  (4.4  through  4.7)  relate  how  the  normal  acceleration 
(a)  should  vary  as  a function  of  the  range  (R),  the  angular  error  (a)  relative 
to  the  desired  straight  line  approach  path  of  slope  (:p)>  and  the  approach 
ground  speed  (-R).  The  set  of  equations  (5.1  through  5.6)  relate  how  the  normal 
acceleration  ('a)  varies  as  a function  of  the  aircraft  input  state  (aT/hi, 

'W^'.  The  function  of  the  autopilot  is  to  manipulate  AT  and  to  cause  Aa 
to  follow  the  desired  command  (■.ar).in  the  presence  of  disturbances  (TW^)  and 
at  the  desired  soeed  of  approach  (-R).  The  landing  technique  most  familiar  to 
all  is  to  manipulate  the  elevens  (A  c)  to  control  descent  attitude  and  hence 
the  soeed  of  approach,  and  to  maniputate  the  throttle  to  control  rate  of  clin.b 
(■,h).  While  this  technique  is  safe  in  that  it  avoids  attempting  to  climb  at 
an  excessive  rate  (e.g.,  no  worse  than  a full  throttle,  constant  airspeed 
climb),  it  is  usua  1 ly unacceptable  for  precise  approaches  due  to  the  sluggish- 
ness of  the  response  of  climb  rate  to  a throttle  input.  This  sluggishness  im- 
pacts both  the  control  response  (e.g.,  how  fast  can  a commanded  ascent/descent 
rate  be  established?)  as  well  as  the  gust  response  (e.g.,  how  fast  can  a gust- 
induced  error  rate  be  arrested?).  The  more  responsive  technique  is  to  manipu- 
late the  throttle  (AT)  to  control  airspeed  and  to  manipulate  the  elevens  (Af^) 
to  control  rate  of  climb  or  more  precisely  normal  acceleration.  One  must  be 
cautious  in  using  this  technique  as  it  is  possible  to  transiently  establish  a 
climb  rate  which  cannot  be  sustained  due  to  insufficient  power  available.  Ad- 
ditionally, if  the  aircraft  is  trimmed  on  the  backside  of  the  drag  curve,  the 
maneuver  is  not  very  stable  (i.e.,  an  increase  in  angle  of  attack  increases 
the  drag  more  than  the  lift.) 


The  behavior  of  normal  acceleration  control  by  elevon  manipulation  will 
be  investigated  by  consideration  of  an  approximate  set  of  differential  equat- 
ions for  which  Au  is  assumed  zero.  A physical  view  of  this  assumption  is  that 
’■‘'o throttle  is  being  deftly  manipulated  to  hold  Au  identically  equal  to  zero 
regardless  of  variations  in  , 'w,  and  Aa. 


From  5. 1 through  5.6, 

“ s -Uo  s 

S(S  -M-) 

Zw  0 


1 

— 

AW 

0 

T 

A0 

Ms 

1 

Aa 

_ _ 

A6^ 


(5.13) 


One  may  compute  the  transfer  function  relating  normal  acceleration  to  elevon 
deflection  as, 


^ ^ -Z^  ^ Z.  Ml  S ^ Up  (Z,  ZJ 

(S  - n*)  (S  - 7J  - Up  Mw 


(5.14) 


Using  the  stability  derivatives  (5.12)  one  nay  assume, 

Zx,  Mv,  <"  Mr  Zw 


(5.15) 
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and  also 


— - - Z,5  + Z(5  M0  S - Zw  Up 

(S  - Mo)  (S  - Z„)  - Uq  Mw 


The  denominator  is  the  most  fundamental  approximation  to  the  short  period  mode 
of  aircraft  motion.  The  steady  state  relationship  between  Aa  and  is  ap- 
proximately, 

— - (5.17) 

Vv 

As  indicated,  one  of  the  main  functions  of  the  autopilot  is  to  reduce  the 
variability  in  how  much  acceleration  is  produced  by  a unit  of  command  input. 
(NOTE:  Only  for  the  case  of  no  feedback  does  one  unit  of  command  input  cor- 

respond to  one  unit  of  elevon  deflection.)  The  above  relationship  (5.17)  im- 
plicitly indicates  the  primary  source  of  this  variability,  being  the  change  in 
My;  produced  by  variations  in  the  location  of  the  aircraft  center  of  mass.  Fig- 
ure 14  depicts  the  behavior  of  the  approximate  equations  (5.13)  in  response  to 
an  elevon  input,  for  two  values  of  M^.  The  use  of  feedback  to  reduce  this  vari- 
ability is  fundamental.  The  quantity  to  be  used  as  tfie  feedback  signal  is,  of 
course,  the  normal  acceleration  itself. 

In  addition  to  reducing  the  variability  of  the  normal  acceleration  re- 
sponse, the  use  of  an  acceleration  feedback  causes  a more  oscillatory  short 
period  mode.  Assume  for  now  that  one  may  neglect  the  effects  of  elevon  lift 
(Z,0  in  (5.13)  and, 

.5^  - {AQq  - Aa ) (5.18) 

Combining  (5.13)  and  (5.18)  one  finds, 

-U  M.  Z K (r 

a _ 0 6 w a 

Aac  ■ (S  - M^)  (S  - Zw)  - Uo  + M^ 

Examination  of  the  denominator  of  (5.19)  shows  that  the  Kg  feedback  has  the 
same  effect  on  the  short  period  mode  as  does  the  "aerodynamic  spring"  or  static 
stability  derivative  (My^).  Increases  in  the  gain  of  the  normal  acceleration 
feedback  will  eventually  have  to  be  dampened  by  a pitch  rate  feedback.  The 
appearance  of  the  term  in  the  numerator  is  due  to  the  acceleration  error 
(Aa^-  - Aa)  formulation  of  (5.18).  This  is,  of  course,  fundamental  to  achieving 
the  desired  relatively  invarient  acceleration  response.  In  the  steady  state, 

,'.a  _ - Up  MiS  Zw  Ka - (5.20) 

Aac  M 0 Zw  - Uq  FV/  - Uq  M(5  Zyj  Kg 


For  large  K^,  variations  in  have  the  desired  negligible  effect  on  the  Aa 
response  per  unit  of  command  input  (Aac).  Using  Ka  = .007  radians  per  feet 
per  second^,  expression  (5.20)  becomes. 


^ = 36 

Aac  1.8+8.5+36 


(5.21) 
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Short  period  response  (At'  = 0)  to  a step 
eleven  input 


Fiqure  15  demonstrates  (1)  the  relative  invariance  of  the  Aa  response  for  two 
values  of  Myy(2),  the  roughly  20  percent  difference  between  the  commanded  ac- 
celeration (Aa^)  and  the  normal  acceleration  (Aa),  and  (3)  the  more  oscillatory 
short  period  mode  caused  by  the  acceleration  feedback  increasing  the  short  per- 
iod frequency  of  oscillation  (Wsp  = 6.8  rad/sec  versus  the  unaugmented  Wsp  = 

3.2  rad/sec).  The  value  of  Aac  in  Figure  15  was  chosen  to  produce  the  same 
initial  elevon  deflection  as  in  Figure  14  (A6r  = 1 degree).  The  decreased 
steady  state  acceleration  achieved  in  Figure  15  is  the  result  of  the  feedback 
signal  dynamically  decreasing  the  elevon  input  from  its  initial  value. 

The  feedback  actually  used  in  Figure  15  was  not  the  acceleration  of  the 
center  of  mass  (Aa), but  rather  the  acceleration  (Aa^)  of  a point  Ij  feet  in 
front  of  the  center  of  mass.  The  initial  jump  in  Aa  is  due  to  a 
lift  upon  the  elevens  (Zg),  neglected  in  the  above  treatment.  The  sense  of 
this  force  is  opposite  to  the  eventual  steady  state  and  is  hence  detrimental 
to  the  stability  of  a Aa  to  Af.^  control  system.  Location  of  the  sensing  ac- 
celerometer at  the  instantaneous  center  of  rotation  (£t  = M6/Z6  = 2.5  feet) 
negates  the  detrimental  effect.  Note  in  Figure  15  that  the  feedback  signal 
'a,D  is  initially  zero  for  £j  = 2.5  feet.  For  the  same  command  input  (-^A^), 
Figure  16  shows  the  amplification  of  the  elevon  lift  effect  by  feedback  of  the 
acceleration  of  the  center  of  mass  (Aap^  = Aa).  For  sufficient  gain  (Kg),  the 
nonminumum  phase  behavior  of  this  feedback  (Aa^,  = .Aa)  would  lead  to  a static 
instability  of  the  short  period  mode.  For  all  remaining  simulation  traces  in 
this  report,  the  sensing  accelerometer  is  located  at  the  instantaneous  center 
of  rotation. 


Figure  17  depicts  the  approximate. short  period  response  with  both  a normal 
acceleration  (Aa^)  and  a pitch  rate  (ag)  feedback.  Using, 


S 

-Mw 

^w 


-Uo  S 

S (S  - M^) 
0 


AW 

r “1 

0 

AG 

= 

^"m 

0 

A6. 


(5.22) 


and, 

A6^  = Kg  (Aa^  - Aa^)  - Kq  A9  (5.23) 

one  finds. 


- Uo  McS  Zw  Ka 

= (S  - Mq  + M^)  (S  - Zw)  -Uq  (Mv,  + Ka  Z„) 


(5.24) 


Not  only  does  the  pitch  rate  feedback  act  to  dampen  the  more  oscillatory  short 
period  mode,  but  it  also  increases  the  steady  state  error  between  Aa^^  and  Aa^-. 
Consider  the  steady  state  of  5.24 


Aa 


m = 


- Up  <a 


Aa,  z„  (Me  ■-  K- 


(5.25) 
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Using  the  stability  derivatives  (5.12)  and  the  selected  feedback  gains  (Ka  = 
- .007  rad/ft/sec^  and  Kg  = - .35  rad/rad/sec)  the  numerical  factors  of  the 
steady  state  relationship  are. 


Aapi  36. 

■ 1.8  + 21.  + 8.5  + 36 

^c 


(5.26) 


For  pure  acceleration  feedback  (‘ap^/Aaf-  = .78)  and  for  pitch  rate  dampened 
acceleration  feedback,  Aa^/Aac  = .53.  Physically,  the  additional  error  is 
developed  because  a prescribed  steady  state  relationship  exists  between  normal 
acceleration  and  pitch  rate  ('.a  = Uq  ■.'").  Not  only  does  the  pitch  rate  con- 
tribute a transient  dampening  effect  but  it  also  introduces  predictable  steady 
state  effects.  This  can  be  countered  by  a feed-forward  path  in  the  autopilot 
(Figure  18)  from  which, 

■'•''e  = ('^c  - + Kq  (AO^,  - AO)  (5.27) 


where 


Equations  (5.27)  and  (5.28)  may  be  combined  to  obtain. 


(5.28) 


(5.29) 


The  feed-forward  path  does  not  affect  system  stability  and  is  equivalent  to 
multiplying  the  input  command  Aac  by  a factor  (1  + Kp/r,a  Uo).  For  the  chosen 
gains  and  approach  speed,  the  numerical  value  of  the  factor  is  1.59.  The  over- 
all idealized  relationship  between  a^^  and  for  the  pitch-rate-dampened, 
acceleration  feedback  would  become  (Aam/Aa^  = .85).  The  feed-forward  feature 
was  not  included  in  this  simulation.  This  is  consistent  with  the  intent  to 
maintain  a conservative  framework  for  this  effort.  If  the  performance  of  the 
coup’’  : guidance/control  dynamics  are  satisfactory  with  Aam'^^.SAaf- then  the  de- 
sign V oncept  is  relatively  sound. 


The  next  autopilot  objective  to  be  addressed  is  the  wind  gust  (AWg)  re- 
sponse of  the  stabilized  airframe.  For  the  landing  mission,  the  gust  charac- 
teristics of  interest  is  the  change  in  flight  path  angle  (y)  produced  by  the 
angle  of  attack  gust  (“g  = AWp/Uo).  indicated  in  section  3,  if  attitude 
is  maintained  (AO  = 0)  during  the  gust,  the  change  in  A.  is  precisely  Artg. 
For  a statically  stable  (M„  negative),  unauqmented  airframe,  the  natural  ten- 
dency is  to  rotate  into  the  wind  and  drive  the  steady  state  angle  of  attack 
(■  -A  = AWa/Uq)  to  zero.  That  is,  in  the  steady  state. 


■ h = Ay  - Aop 


(5.30) 


which  in  the  steady  state  equals, 

Afl  _ '-'o 

■,  .G  ■ Zw  M-  - Uo  Mw 


(5.33) 


(5.34) 


For  the  stability  derivatives  of  5.12,  then  AO/Aa^  = .9  in  the  steady  state. 
The  response  of  the  basic  airframe  to  a downqust  (\-<q  = - 3.4  degrees  or 
'W'"  = - 5 ft/sec)  is  shown  in  Figure  19a  and  confirms  expression  (5.34).  As 
seen  earlier,  normal  acceleration  and  pitch  rate  feedbacks  to  the  elevens  may 
be  treated  as  increments  to  and  M;'j,  respectively. 


^w  ^ ^a  ^ i ^w 


(5.35) 


The  conclusions  are  quite  clear.  The  acceleration  feedback  increases  the  ten- 
dency to  rotate  into  the  wind  and  is  beneficial  to  the  gust  response,  while 
the  pitch  rate  feedback  inhibits  rotation  into  the  wind  and  is  detrimental  to 
the  gust  response.  Figure  19(b)  and  Figure  20(a)  demonstrate  these  individual 
characteristics.  In  keeping  with  the  aforementioned  conservative  design  ap- 
proach, the  chosen  values  for  Ka  and  make  the  combined  augmented  system 
somewhat  more  gust  sensitive  in  the  steady  state  (Figure  2Cb)  than  the  basic 
airframe  (Figure  19a),  but  still  substantially  less  than  the  A9  = 0 analysis 
of  section  3. 


The  autopilot  discussion  thus  far  has  been  concerned  with  control  of  the 
short  period  or  anale  of  attack  motions  of  the  airframe.  The  discussion  and 
simulation  results  have  proceeded  on  the  assumption  that  thrust  variations  of 
the  engine  maintain  AU  identically  equal  to  zero.  In  a similar  fashion,  phugoid 
airframe  motions  can  be  visualized  by  the  assumption  of  AW  = 0.  This  is  equiv- 
alent to  the  airframe  having  a quick  stable  short  period  mode  so  that  angle  of 
attack  motions  may  be  neglected.  The  phugoid  mode  is  dominated  by  variations 
in  the  X and  Z forces  due  to  AU  and  AO  perturbations. 
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Figure  19.  Short  period  gust  response  (AWq  = -0.5  ft/sec) 
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Figure  20.  Short  period  gust  response  (AWq  = =5  ft/sec) 
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(5.37) 


The  liqhtly  dampened  approximation  to  the  phugoid  mode  is  shown  in  Figure  21. 

The  \l  = 0 constraint  was  approximated  on  the  simulation  trace  by  increasing 
(-  1.  versus  -.1).  The  most  direct  method  of  damping  the  phugoid  is  to 
increase  the  drag  of  the  airframe.  Modulation  of  the  thrust  of  the  engine  is 
equivalent  but  too  small  in  magnitude  to  be  useful.  Deployment  of  drag  brakes 
would  permit  steeper  descents  as  well  as  dampen  the  phugoid  but  this  approach 
is  not  considered  here.  The  usual  means  of  increasing  phugoid  damping  is  to 
^eed  low  frequency  oitch  attitude  information  back  to  the  elevens.  (Gust 
response  considerations  preclude  high  frequency  pitch  attitude  feedback).  Be- 
cause the  RPV  mission  tends  to  require  extended  periods  of  circling,  a verti- 
cal gyro  is  not  usually  part  of  the  RPV  sensor  package  and  low  frequency  filter- 
ing of  pitch  rate  to  obtain  quasi-attitude  information  is  undesirable  because 
of  the  filter  initialization  difficulties  which  are  typical  of  such  implemen- 
tations. The  approach  taken  here,  therefore,  is  to  rely  on  the  Aa,  A6  auto- 
pilot feedbacks  as  well  as  the  low  frequency  A,  X guidance  feedbacks  to  main- 
tain control  of  the  phugoid  motions. 


Figure  22  is  a simulation  trace  of  the  complete  airframe  motions  (5.1  to 
5.12)  in  response  to  a 5 ft/sec  wind  gust  (AWg).  The  phugoid  mode  has  less 
damping  than  the  classical  phugoid  (AVJ  = 0)  due  to  coupling  of  the  AW  motions 
(”w  = - .1)  with  the  AU,  A6  motions.  Theperiodqf  the  phugoid  motion  (Wp”^  = 

13  seconds  is  relatively  unaffected  by  the  AW  degree  of  freedom.  The  in- 
corporation of  a normal  acceleration  feedback  affects  both  the  short  period 
and  the  phugoid  modes  of  motion.  As  indicated  in  f-igure  23,  the  damping  of 
the  phugoid  mode  is  more  like  the  classical  phugoid  case  (Figure  21).  The  ac- 
celeration feedback  has  quickened  the  short  period  mode  and  made  the  aW  = 0 
approximation  of  Figure  21  more  valid  than  for  the  unaugmented  airframe  case 
of  Figure  22.  Another  effect  of  the  feedback  is  to  convert  the  phugoid  motions 
into  lower  frequency  (Wp"'  = 28  seconds),  slightly  higher  amplitude  motions. 

This  is  beneficial  as  the  lower  frequency  behavior  is  more  amenable  to  stabil- 

ization by  the  guidance  feedbacks  (X,X)  on  approach.  A pure  pitch  rate  feed- 
back has  the  same  low  frequency  effects  upon  phugoid  stability  as  does  the 
normal  acceleration  feedback  (e.g.,  Aa  ==  Po  important  difference, 

however,  is  the  magnitude  of  the  excitation  of  the  phugoid  motion.  As  can  be 
seen  in  a comparison  of  Figures  23  and  24,  the  A9  induced  by  the  initial  wind 

gust  is  much  greater  for  the  Aa  feedback  than  for  A0  case.  It  is  apparent 

that  the  price  to  be  paid  for  the  AWg  gust  alleviation  character  of  the  Aa 
feedback,  is  increased  excitation  of  the  phugoid  mode.  The  combined  effects 
of  the  A0  and  Aa  feedbacks  are  shown  in  Figure  25. 

The  coupled  set  of  perturbation  equations  (5.1  through  5.6)  together  with 
the  stability  derivatives  (5.12)  are  stable.  As  a result  of  being  stable, 
when  these  equations  are  perturbed  from  their  trim  value  by  a step  wind  gust 
(AWg),  the  aircraft  will  eventually  return  to  this  trimmed  aerodynamic  state 
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Figure  22.  Coupled  'W 
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Figure  23.  Augmented  response  (K, 
5 ft/sec) 
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(AV  = ‘W-  = ' • = 0).  'Ihis  behavior  is  clearly  shown  in  Figure  22  where  one 
Should  r^vall  that  ' = ' 'J  + 'Wg  and  that  being  nonzero  is  consistent  with 
the  stated  aerodynamic  trim  of  the  aircraft.  The  incorporation  of  an  autopilot 
mav  or  may  not  modify  the  behavior  described  above.  The  autopilot  feedback 
(5.2.T)  is  transient  in  nature  so  that  the  stabilized  airframe  (Figure  25)  also 
returns  to  the  original  aerodynamic  trim.  If  the  autopilot  incorporates  a 
pitch  trim  integrator  such  as, 

t 

' t “ ^a  ( '^c  ■ '^)  ■ /p  (■'•^c  ■ dt  (5.3R) 

then  the  long  term  behavior  of  the  aircraft  is  modified  as  is  nonzero  in 
the  steady  state  (Figure  26).  It  is  the  opinion  of  this  author  that  caution 
be  exercised  in  the  use  of  low  freouency  feedbacks  in  a necessarily  responsive 
autopilot  channel.  Crror  or  trim  integrators  such  as  5.38  are  the  least  offen- 
sive form  of  rhese  type  feedbacks  but  the  general  feeling  still  applies.  The 
feedbacks  can  have  a beneficial  effect  on  the  long  term  stability  of  the  air- 
frame (e.q..  note  the  increased  phugoid  damping  of  Figure  26  compared  with 
rigu*'e  2-'),  but  such  effects  can  also  be  obtained  by  guidance  phenomena  (', 
etc.).  To  retain  the  conservative  framework,  expression  5.38  is  used  as 
the  acceleration  autopilot  feedback  for  the  remainder  of  the  simulation  traces 
(see  Appendix  A for  the  complete  autopilot  representation).  Figure  27  de- 
picts the  gust  response  of  the  complete  autopi lot/airframe  system. 

6.  COUPLED  GUIDANCE  AND  CONTROL  FUNCTIONS  ' 


The  approach  taken  thus  far  has  been  to  treat  the  guidance  ^unction  and 
the  control  function  in  an  independent  fashion.  This  does  not  imply  that  the 
potential  of  coupling  between  the  guidance  dynamics  and  the  control  dynamics 
has  been  ignored.  The  near-field  gain  limiting  was  incorporated  into  the 
guidance  law  to  achieve  compatibility  between  the  guidance  dynamics  and  the 
autopi lot/ai rfname  dynamics.  The  simulation  effort  will  be  completed  by  the 
merging  of  these  two  functions. 

Figure  28  is  similar  to  Figure  13  with  the  exception  that  the  commanded 
normal  acceleration  (‘ac)  and  the  commanded  airspeed  ("Uc)  are  filtered  by  the 
airframe./autopi lot  dynamics  before  computation  of  the  guidance  kinematics  (■, 
i).  Since  the  guidance  gain  is  not  limited  (R^  = 0),  the  expected  near-field 
flight  path  instability  is  produced.  For  the  autopilot  implementation  of  Fig- 
ure A-4,  the  range  at  which  the  instability  occurs  is  on  the  order  of  300  feet 
for  the  chosen  n = 15.  One  might  now  expect  to  see  a stability  analysis  of  the 
guidance  law  (4.4)  and  the  U = 0 approximation  to  the  stabilized  airframe 
(5.19).  Unfortunately , this  stability  analysis  is  not  so  easy.  The  system  of 
equations  depicted  by  Figure  28  is  linear  , but  due  to  the  time-varying  coef- 
ficients produced  by  the  constantly  decreasing  range,  an  analytical  treatment 
is  really  quite  difficult  and  fortunately  unnecessary.  The  near-field  gain 
control  eliminates  the^analytic  difficulties  by  conversion  of  the  measured  an- 
gular information  (\,  ) to  its  equivalent  (hf.  He)  and  stabilization  of  the 

system  gain  at  the  R^i  value.  More  importantly,  proper  adjustment  of  the  gain 
control  will  be  seen  to  eliminate  the  instability  along  with  the  associated 
analytic  difficulties. 
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Figure  ?9.  Coupled  control  response  ( 
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Figures  30  and  31  depict  the  gust  response  of  the  coupled  system  with  gain 
limiting  (R|.,  = 300  ft).  As  in  Section  3,  the  system  is  perturbed  (Rq  = 600  H) 
from  its  initial  trim  on  the  glide  slope.  The  instability  is  still  in  evidence 
though  its  severity  is  bounded  by  the  gain  limiting  (R^i  = 300  feet  to  touchdown. 
To  provide  some  gain  margin  for  the  autopilot  configuration  of  Figure  A-4,  the 
design  value  for  the  gain  parameter  (R„,)  was  chosen  to  be  1,000  feet. 

The  coupled  response  of  the  design  conf iguration  (n  = 15,  Rm  = 1,000  ft) 
is  given  in  Figures  32  and  33.  The  idealized  guidance  response  (,'.a  = ,'ac)  was 
previously  given  in  Figure  13.  The  most  apparent  difference  between  the  two 
figures  is  the  low  frequency  drift  (Aa(~,  Aa^)  of  Figure  32.  This  phenomenon 
is  the  result  of  the  trim  integrators  previously  mentioned  and  more  quantita- 
tively de()icted  in  the  ' '[t  trace  of  Figure  33.  Comparison  of  the  .‘.a  and  ,'.3^ 

traces  demonstrate  that  a successful  approach  is  not  dependent  upon  holding 
ap  near  zero.  As  a matter  of  interest,  the  reader  should  also  note  the  effec- 
p've  control  of  descending  flight  path  angle  (jy|  -s:  8°)  by  command  limiting  of 


Figures  34  through  37  depict  the  coupled  guidance  and  control  responses  of 
the  design  configuration  to  a AWg  gust  input  (Rn  = 600  and  300  feet,  respec- 
tively. In  both  cases,  the  peak  "hf  was  one  (1)  foot.  The  invarient  Ah£  per- 
turbation is  consistent  with  the  concept  of  a constant  gain  system  for  all 
ranges  less  than  R^,  = 1,000  feet.  The  worst  case  computations  of  Section  3 
( = 0)  indicated  a peak  ',hE  of  1.5  feet  (e.g.  , See  Figures  11  and  12).  The 

difference  is  due  to  the  gust  alleviation  properties  of  the  stabilized  airframe. 

7.  SUMMARY  OF  CONCLUS IONS 


The  efforts  of  the  preceding  six  sections  have  tried  to  indicate  a sug- 
gested direction  for  the  design  of  an  RPV  approach  guidance  and  control  system. 
Every  effort  was  made  to  maintain  a conservative  framework  for  the  nominal  de- 
sign configuration  developed  within  those  six  sections.  As  a result,  it  is 
felt  that  the  performance  of  the  concept  can  be  significantly  improved  by  a de- 
tailed stability  and  control  analysis,  particularly  in  the  near-field  region 
with  the  intent  of  selecting  optimum  values  for  the  parameters  (Rm,  n,  Ka , K§) 
for  use  with  a specific  RPV  airframe/engine  design.  The  more  pertinent  results 
of  the  eftort  of  the  first  six  sections  follow. 

a.  An  approach  guidance  law  was  developed  which  was  felt  to  have  applica- 
tion to  the  general  problem  of  IFR  landing  of  aircraft  of  any  type.  For  contin- 
uous measurement  of  the  angular  deviation  of  the  aircraft  from  the  glide  path 
(:),  its  time  rate  of  change  [a]  and  the  horizontal  range  (R)  of  the  aircraft 
from  the  landing  site,  the  guidance  law  was  shown  to  have  inherent  stability 
even  as  the  range  approaches  zero.  To  achieve  this  inherent  stability,  the  band- 
width of  the  guidance  law  is  reguired  to  vary  as  R"^ . 

h.  To  compensate  for  the  effects  of  lags  in  the  response  of  the  airframe 
to  the  guidance  law  demands,  a near-field  gain  control  was  incorporated.  This 
gain  control  was  shown  to  limit  the  bandwidth  of  the  guidance  law  for  R<  Rm  to 
a value  compatible  with  the  airframe's  ability  to  respond.  The  guidance  law 
bandwidth  must  also  be  compatible  with  the  data  rate  and  noise  spectra  associa- 
ted with  measurement  of  A,  a,  and  R.  In  many  applications,  a direct  measurement 
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Figure  32.  Coupled  guidance  response  (n  = 15,  R^i  = 1000  ft 
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Figure  33.  Coupled  control  response  (n 
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Figure  35.  Coupled  control  response  (n  = 15,  Rm  = 1000  ft 
= 40,  ",Wr  = 5 ft/sec,  Rd  = 600  ft) 


Fiqure  37.  Coupled  control  response  (n 
,p  = 40,  aWq  = 5 ft/ sec,  R[) 
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of  will  not  be  available.  System  performance  (i.e.,  guidance  law  bandwidth) 
will  ^h^refore  be  influenced  by  one's  ability  to  approximate  ■ within  the  de- 
sired frequency  region  by  operating  upon  the  measured  value  of  A. 

c.  The  near-field  (R  < R^)  gust  response  characteristics  of  the  guidance 
law  were  quantitatively  related  to  the  guidance  parameters  R^  and  n,  and  equiva- 
lently, to  the  bandwidth  of  the  guidance  law.  This  analytic  treatment,  which 
neglected  airframe  dynamics,  was  shown  to  represent  a worst  case  situation  due 
*•0  the  gust  alleviation  properties  of  the  stabilized  airframe. 

d.  Command  limiting  (A^)  within  the  guidance  law  was  shown  to  provide  an 
effective  control  over  the  maximum  angle  of  descent  when  a measurement  of  flight 
path  angle  (. ) is  not  available.  The  command  limiting  is  necessary  to  avoid  a 
build-up  of  airspeed  under  flight-idle  engine  conditions.  The  technique  may 
also  be  applied  to  climbing  flight  to  avoid  conditions  of  insufficient  power 
available,  but  this  application  was  not  considered  in  this  simulation  effort. 

-.  A very  simple  "approach  autopilot"  concept  was  suggested  for  the  mini- 
RPV  landing  function.  The  autopilot  relies  on  normal  acceleration  and  pitching 
rate  feedbacks  to  control  and  stabilize  the  short  period  and  phugoid  modes  of 
airframe  motion. 


An  attempt  has  been  made  within  this  document  to  carefully  distinguish 
between  the  guidance  function  and  the  control  function.  The  intent  was  to 
clearly  discriminate  the  design  objectives  of  each.  When  one  considers  the 
physical  mechanization  of  the  two  functions,  however,  there  is  no  associated 
implication  that  the  two  be  clearly  discernable  (e.g.,  guidance  calculations  in 
a ground-based  computer,  control  calculations  within  the  air  vehicle).  Most 
certainly,  the  minimum  required  airborne  computation  is  that  associated  with 
the  control  function.  With  the  advent  of  microprocessor  technology,  however, 
one  is  inexorably  drawn  toward  centralization  of  the  guidance  and  control  calcu- 
lations in  what  should  be  referred  to  as  an  "Approach  Guidance  and  Control  Com- 
puter." The  inputs  to  the  computer  would  be  communicated  A and  R from  the 
ground  based  radar  or  optical  tracker,  and  the  airframe  normal  acceleration  and 
pitching  rate  from  the  airborne  sensors.  The  outputs  of  the  computer  would  be 
electrical  signals  to  drive  the  control  surface  sensors  (i.e.,  a "digital  auto- 
pilot") for  an  RPV  application,  or  to  drive  a flight  director  for  the  IFR, 
piloted  application. 


APPENDIX 


This  appendix 
frame,  auto-pilot. 


contains  the  scaled  analog  simulation  diagrams  for  the  air- 
and  associated  kinematic  relationships. 
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